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Abstract 

We present a gauge invariant formalism to study the evolution of curvature perturbations during 
decay of the inflaton or some other field that dominates the energy density. We specialize to the 
case where the total curvature perturbation arises predominantly from a spatially fluctuating decay 
rate. Gaussian fluctuations in the decay rate are present if the inflaton coupling or mass depends 
on some light field, which fluctuates freely during inflation. The isocurvature fluctuations of the 
light field seeds the total curvature perturbations, whose amplitude freezes after inflaton decay. We 
also analyze the resulting curvature perturbation starting with non-negligible values of the initial 
inflaton curvature. Lastly, we study the effects of the energy density stored in the light field. 



I. INTRODUCTION 



Inflation is the main contender for an explanation of the observed adiabatic density 
perturbations with a nearly scale invariant spectrum [lj]. However, recently, alternative 
mechanisms for generating the density perturbations have also been discussed. In the cur- 
vaton scenario, isocurvature perturbations of some light "curvaton" field are converted into 
adiabatic perturbations in the post- inflationary universe Q]. 

Another interesting proposal is that the perturbations could be generated from the fluc- 
tuations of the inflaton coupling to the Standard Model degrees of freedom HO, It has 
been argued that the coupling strength of the inflaton to ordinary matter or the inflaton 
mass, instead of being a constant, could depend on the vacuum expectation values (VEV) of 
the various fields in the theory. If these fields are light during inflation their quantum fluc- 
tuations will lead to spatial fluctuations in the inflaton decay rate. As a consequence, when 
the inflaton decays, adiabatic density perturbations will be created because fluctuations in 
the decay rate translate into fluctuations in the reheating temperature. A variation of the 
above scenario is that the inflaton decays into heavy particle, whose decay is governed by a 
varying decay rate. If these particles come (close) to dominating the energy density of the 
universe, curvature perturbations are produced The light fields whose fluctuations leads 
to a fluctuating decay rate could be normal matter fields, such as the flat directions in the 
minimal supersymmetric Standard Model (MSSM), or a right handed neutrino []|. 

In this paper we will discuss the evolution of the curvature perturbation during decay, 
in the case of a fluctuating decay rate. We will use the gauge-invariant formalism for 
metric perturbations as developed by Bardeen for a review see Ref. The curvature 
perturbation on large scales (defined on a uniform density spatial hypersurface) ( remains 
constant for purely adiabatic perturbations. But ( can change on large scales due to a non- 
adiabatic pressure perturbation. In a multi-fluid system such pressure perturbations can be 
important in the presence of relative entropy perturbations. Said in another way, entropy 
(or isocurvature) perturbations can feed the curvature perturbation. Therefore, to compute 
the final perturbation spectrum in a multi fluid system it is necessary to keep track of the 
evolution and perturbations of various fluids. 

In the case of a fluctuating decay rate, the multi fluid system consist of the decaying field 
(the inflaton, or some other field which has dominant energy density), the radiation bath, 
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and the light or flat direction field that is responsible for the fluctuating decay rate. The 
relative curvature perturbation between the inflaton and the flat direction acts as a source 
for the total curvature perturbations. This can be understood intuitively from the fact that 
fluctuations in the inflaton decay rate leads to fluctuations in the reheat temperature of 



the universe, given by T rh ~ X^jm^Mp, where is the mass of the inflaton and Mp = 
2.436 x 10 18 GeV is the reduced Planck mass. The fluctuations in T can be translated into 
fluctuations in the energy density of a thermal bath with 5p 7 /p 7 = — (2/3)<5r/r [3]. The 
factor 2/3 appears due to red-shift of the modes during the decay of the inflaton whose 
energy still dominates. 

The inflaton decay rate is T ~ m^X 2 . The decay rate fluctuates if either A or is a 
function of a fluctuating light field. The former case can e.g. arise in the MSSM through 
terms in the superpotential of the form 

WBX^HH + ^qq + ^q.q. + ^—qq., (1.1) 

with the inflaton field a standard model singlet, H and H the two Higgs doublets, and 
q and q c the quark and lepton superfields and their anti-particles. M is some cut-off scale 
which could be the GUT scale or the Planck scale. In this paper we will not be concerned 
with the particulars of the decay products and the light fields. What we will take from the 
above example though, is that there can be two kinds of couplings: 

A I A o (l + f + •••) , direct decay. 
I A ° (l? + •■■) ' indirect decay, 

with the ellipses standing for higher order terms. In the MSSM example, the "direct" decay 
mode corresponds to inflaton decay into Higgs fields, whereas the "indirect" decay mode 
corresponds to decay into (s) quarks and anti-(s) quarks. S is the expectation value of a 
light fluctuating field. Effective couplings of this form can result from integrating out heavy 
particles. Another possibility is that the inflaton, or some other field that dominates the 
energy density of the universe (e.g. non-relativistic inflaton decay products), has an effective 
mass generated through a coupling to the flat direction: = XS. The fluctuation in the 
decay rate for the various cases is 



6T 

r 



2|f, direct decay. 

2^, indirect decay. (1.3) 
fluctuating mass. 
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II. BACKGROUND EQUATIONS 



The background equations of motion governing the dynamics of the inflaton, the radiation 
bath, and the flat direction field in a Robertson Walker universe are 

p = - p<j> (3H + T), (2.1) 

p 7 = -AH Pl + T P(t>) (2.2) 

p s = -3H Ps (l + u). (2.3) 

where a dot denotes differentiation w.r.t. coordinate time. The Hubble parameter H = a/a, 
with a the scale factor of the universe, is given by 

# 2 = ^r(P^ + P7 + Ps)- (2-4) 
Here u = ps/ ps is the equation of state parameter for the flat direction field S. When 
H 3> wis, the flat direction field is over damped and remains effectively frozen, and uj — ► — 1. 
For H < ms, the flat direction field oscillates in its potential. For a quadratic potential 
lu = averaged over one oscillation. To determine u and the decay width T(S) we need the 
evolution of S: 

S + 3HS + V S = (2.5) 

with V s = OV/dS. Then 

Ps S 2 /2 + V(S) 
It is useful to introduce the following dimensionless parameters 

il a = —, s = —, g = , h = . (2.7) 

p m 

We leave the mass scale m unspecified for now, it can be taken m 5 . Note however that if 
finite temperature effects or effects of non-renormalizable operators are taken into account, 
then m s is time-dependent and we cannot set m = (ms) e s- Further, we eliminate f2 7 from 
all equations by use of the Friedman constraint: 

n s + + fi 7 = i (2.8) 

The dimensionless background equations are then 

% = + ^(3cu - 1)) (2.9) 

Q' s = n s ((l - 3cu) - % + (3cu - l)n s ) (2.10) 

ti = |( - 4 + % + (1 - 3cu)n s ) . (2.11) 



and 

' + (H' + W- a (2i2) 

Here, prime denotes differentiation with respect to the number of e-foldings N = In a. 

III. GAUGE INVARIANT PERTURBATIONS 

We are interested in the long wavelength regime of the perturbations where the comoving 
scale is much larger than the Hubble horizon. Following Ref. s| we define the curvature 
perturbation ( on spatial slices of uniform density p with the line element 1 

ds 2 = a 2 (t) (1 + 2C) 5 ij dx i dx j , (3.1) 

where a is the scale factor. The time evolution of the curvature perturbation on large scales 
is [lQ, 11, Q 

C = ~jfp8P*»d, (3-2) 

where P nac j = SP — c 2 Sp is the non-adiabatic pressure perturbation. The adiabatic sound 
speed is c 2 s = P / p, where P and p are the total pressure and energy density respectively. 
For a single field 5P na d = 0, and therefore on large scales the curvature perturbation is pure 
adiabatic in nature with ( = constant. 

The total curvature perturbation Eq. ()3.2j) can be written in terms of the various com- 
ponents 

c = E-c- (3-3) 

a P 

In our case a = <fi, S, 7. Isocurvature or entropy perturbations describe the difference be- 
tween the curvature perturbations (n| 

Saf, = HCa~ Cfd = -3# f ^ - ^ V (3.4) 

V Pa Pd J 

With the help of Eqs. (j3.3l3.4J) . we can obtain a useful relationship 

C* = C + \ E > (3-5) 

6 (3 P 



1 Linear perturbations around a flat FRW metric can be denned by the line element ds 2 = — (1 + 2<fi)dt 2 + 
2aB_idtdx l + a 2 [(1 — 2ip)5ij + 2E.ij] dx % dx\ where we have used the notation of Ref. 9] for the gauge 
dependent curvature perturbation ip, the lapse function cj>, and the scalar shear \ = a2 E ~ a B- The 
quantity £ is related to the curvature perturbation ip, on a generic slicing, by £ = —ip — HSp/ p. 
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where a, (3 = </>, S, 7. 

In the presence of more than one field the non-adiabatic pressure perturbation can be 
non-zero <5P na d ^ 0, and therefore the total curvature perturbations evolve in time. The non- 
adiabatic pressure has two contributions, one is from the intrinsic pressure perturbations and 
the other from the relative pressure 

5-Pnad = 3Pintv,a + SP re \ , (3.6) 

a 

where SP intTta = 8P a — c a $Pa, and c 2 a = P a /p a is the adiabatic speed of sound for a respective 
fluid. Note that the intrinsic pressure perturbation vanishes for a fixed equation of state. 
Thus, <5Pi nt r i7 = 0, and (during the stage of inflaton oscillations) <5-Pmtr,</> = 0. For a fixed u 
parameter, also the intrinsic pressure perturbation for the flat direction vanishes <5Pi ntrj s = 0. 
However the relative pressure perturbation is non-zero and can be expressed in terms of the 
various entropy perturbations as ^j] 

QHp af3 v pi 



Using Eqs. ()3. 213. 313. 413.7)) . the evolution of the total curvature perturbation becomes 

h n . . „ i, 

J 2 

During the decay of the inflaton in relativistic species there is non-adiabatic energy trans- 
fer from the inflaton to the radiation bath, see Eqs. ()2.ip2.2|) : 



£ = ~2 (^P4>P7 S <h ~ ( u ~ ^)PsP~/Ss~/ + up^psS^ . (3.8) 



= -Tpd, (3.9) 
Q 7 = Tpf. (3.10) 

Due to fluctuations in T, the non-adiabatic energy transfer is also subject to perturbations. 
Therefore, the evolution of the individual curvature perturbations and £ 7 is determined not 
only by non-adiabatic pressure perturbations but also by the non-adiabatic energy transfer 
perturbations 5Q n ad,a- Analogous to the non-adiabatic pressure perturbations, these can be 



split in an intrinsic and relative part 
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^Qnad,a = ^Qintr.a + $Qrel,a ■ (3.11) 



The intrinsic non-adiabatic energy transfer is defined as 



$Qmtv,a = 5Q a ~ ^5p a , (3.12) 
Pa 



whose contribution is zero if Q a is a function of local energy transfer 5Q a = Q a Sp a /p 
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12 1 . The non-zero contributions are 



r 

8Qmt,<j> = —7^P<i>S<f>s , (3.13) 

f r 

5Qint, 7 = -^jjP^S^S — T^jPcpStfrf , (3-14) 

where we have used <5r/r = 5S/S, see Eq. (|1.2jl . to write the exp ressions in terms of 5*0,5. 
The relative non-adiabatic energy transfer is given by [ill . I12I ] 

5Qrcl,a = J^PpSaP ■ (3.15) 

The non-zero components are 

8Qxe\,(t> = ^^(PSS^ S + PyS^y) , (3.16) 

5Qrel,7 = —^^{P4>S^ + PsS 7 s) . (3.17) 

The evolution of the individual gauge invariant curvature perturbations (on uniform density 
hypersurfaces) is fed by the non-adiabatic perturbations through 

• 6~H 2 5P mtra H5Q na d, a /„ , i 

C, a — : : • {6.18) 

Pa Pa 

The evolution equations of the curvature perturbations then are 

* = & (r-fW-^, (3-19) 




^ = t(- f + ^)^ + ^( 1 -0) S *" (3 ' 20) 

• Oil Ol\ n ir Q , . 

( s = EE£. 3.21 

PS 

The equation for (s needs some explanation. First of all, note that the light field S does 
not undergo energy transfer, unlike the inflaton and the radiation bath. If the equation of 
state parameter is a constant, as it is for example during S oscillations, the non-adiabatic 
pressure perturbation vanishes and (s remains constant. This just reflects the fact that 
SS/S oc Sps/ps is constant for a quadratic potential, as S and 5S follow the same equation 
of motion. 

During slow roll the equation of state parameter u is changing with time and the non- 
adiabatic pressure perturbation is non-zero. If the flat direction field is rolling in a potential 



dominated by non-renormalizable terms, this can lead to appreciable damping of (SS/S), 
as discussed in Ref. p. Therefore, we will assume a quadratic potential for S with non- 
renormalizable terms negligible small. This assumption is justified if the VEV of S is small 
enough, \m 2 s S 2 <C A 2 S 2n ~ 2 / 'M 2n ' 6 , with n > 4. Small VeVs are naturally obtained in low 
scale inflation. 

However, this still prohibits us from taking the 'frozen' limit S — > in a naive way. The 
problem is that in this limit uj — ► —1, and £s — > oo becomes ill defined. To do things 
properly one has to take the slow roll into account. We write uo = — 1 + e, with e = 2T/V in 
the slow roll limit (S — > 0, and T <C V). Here T = ^S 2 the kinetic energy, and V = \m 2 S 2 
is the potential energy. Then 

Pmtr,s = —ops = ep s . (3.22) 

ps Ps 

During slow roll the curvature perturbation of S is slowly changing with time. 
Eliminating £ 7 and p 7 we can rewrite Eqs. (|3.19I3.2(J|) in dimensionless form 

c ' = " -4 + n, + 0,(1 - { ( - "*< 3 + + *W + - D)c 

-3f! s (l + w)(3« - l)Cs + (3 + s//»)!i*C*} , (3-23) 
a =(-4 + !i , + !is( l-3a,)) 5 ^C + ^C S 

+ ?T3A ( " 9 ' + I ( " 4 + + 05(1 " 3lj)) ) C * ' (3 - M) 
where as before a prime denotes differentiation w.r.t. the number of e-foldings N = In a. 



IV. PHYSICAL SITUATIONS 



To simplify and explore the properties of the evolution equations for the perturbations 
we will discuss two special cases in this section. 



A. S non-dynamical 

First we consider the case where the field S is strictly frozen and is non-dynamical. In 
addition we require S to have negligible energy density. The evolution of the perturbation 
is governed by the fluctuating inflaton decay rate, irrespective of what is the cause for 
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FIG. 1: C as a function of N = In a for different initial inflaton perturbations. The bottom 
line corresponds to Q>(0) = (solid). The lines above are for Q>(0) = 0.001 (dashed-dotted) , 
C</>(0) = 0.05 (dotted), and C</>(0) = 0.1 (dashed) respectively. 



the fluctuations. Note that although becomes ill defined in this limit, the combination 
g'Cs — <5r remains finite. The decay rate is independent of time. It is then useful to introduce 
the quantities 



f = -= 9 - 
J H h 



ST 

Or = — ■ 



In the non-dynamical limit the background equations become: 



/' = -{(-4 + 0^). 



(4.1) 

(4.2) 
(4.3) 



The perturbation equations are 

c = 

(<t> = 



(3 + 

-4 + fi 
-4 + fi*)/ 



(C* - 
(C* - 



(4.4) 
(4.5) 



2(3 + /) - 3 + / 

We solved the equations numerically. The results are shown in Fig. 1. Plotted is the 
ratio £ = C/Sr- This quantity is independent of the value of 5r, as expected |3J]. The solid 
line corresponds to negligible initial curvature perturbation of the inflaton: (^(No = 0) = 0, 
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where we have introduced the notation ( a = ( a /5r- In all plots, /(0) = 10 2 , i.e., the 
evolution starts at H = 10 2 T. H ~ T occurs at JV « 3. We find C « -0.166 = -(1/6). 

n 

This agrees with the results of 3], who found for the gravitational potential during radiation 
domination if) = — (2/3)C = (l/9)<5 r 2 • 

Furthermore, we looked at the effects of a non-zero initial inflaton curvature. For C/>(0) ^5 
10 -3 the final curvature is indistinguishable from ($(0) = 0. But for larger initial curvatures 
deviation from ( = — <5r/6 arise. For Cp(0) = 10 -2 the deviation is ~ 5%, whereas for 
C*(0) = 5 x 10~ 2 it is ~ 30%. 



B. S oscillating 

It is possible that the inflaton decays while the flat direction field is oscillating. This is 
the case for T < ms- Averaged over one oscillation u = 0, and the averaged decay rate is 
constant in time. The dimensionless background equations for this case are 

% = % ( i - f - - n s ) , (4.6) 



h 

n' s = n s (i - ^ - n s ) , (4.7) 

f = -£(-4 + 0^ + Q 5 ). (4.8) 

The perturbation equations are 

>, (-^(3 + /) - 3fl s ) C + 3QgG + (3 + /)^C0 r , p x 

C, - (C* - - 3^ , (4-10) 

= 0. (4.11) 

Since 5 has a fixed equation of state in this case, (s is constant, and it is straightforward to 
incorporate the evolution of S. It is easy to see that when Qs and Qs Cs are negligible, the 
above equations are the same as in the non-dynamical limit, and the result ( = — Sr/Q is 
obtained. This is shown by the solid line in Fig. 2. For Q$ large enough, the final curvature 
starts to deviate from this value. A large fls is possible especially in the scenario discussed 
in 0| , where it is not the inflaton but its decay products that have a varying decay rate. 



2 The metric potentials ip and £ are related to each other by C = — (2/3) (i? 1 ip + V0/(1 + Lit) — ip- For 
constant tp and w — 1/3, corresponding to radiation domination, ( — —(3/2)?/). For details see Ref. 
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FIG. 2: C as a function of N = In a for different contributions of Qs an d Cs- The solid line 
corresponds to negligible small, the dotted line corresponds to £ls(N = 0) = 10~ 5 and Cs = 10 2 
(for direct decay), and the dashed line corresponds to f2s(iV = 0) = 5 x 1CP 3 and (s = 1/3 (indirect 
decay). In all plots 6r = 10 -5 . 



The field S does not undergo any transfer of energy: Q$ = 0. As a result, a non- 
negligible energy density Qs can only affect the total curvature perturbation as a result 
of non-adiabatic, relative pressure perturbations, see Eq. (|3.7j) . Since we assume that the 
inflaton red shifts as non-relativistic matter after inflation, there is no relative pressure 
perturbation between the inflaton and the S field: S^s = 0. There is a relative pressure 
perturbation though between the radiation bath and S, and S y s 0- However, it is only 
after inflaton decay that the energy density stored in the radiation bath becomes appreciable, 
and that Qs increases. Then, the relative pressure perturbation leads to the conversion of 
isocurvature perturbations to curvature perturbations. This is the idea behind the curvaton 
scenario jQ]. The final curvature is approximately given by the function of both £ 7 and (s'- 



C w (1 - r)C 7 + r( s 



(4.12) 



with 



4fi 7 + 3^5 



(4.13) 



S decay 



Here r is to be evaluated at the time of S-decay. Since S does not undergo energy exchange 
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it acts like a spectator field during inflaton decay, and £ 7 ~ — <5r/6, independently of the 
value of Qs and decay channels, i.e. direct /indirect decay modes, see Eq. (jl.2D . 

The curvature perturbation of the S'-field, £5, can be much larger for direct decay than for 
either indirect decay or the fluctuating mass case, see Eq. (jl.3|) . The curvature perturbation 
of the flat direction field during the epoch of S'-oscillations is 

Cs = ^ = ^. (4^14) 

dps ob 

For direct decay 5r = ST/T = 25S/M < (5, where the last inequality follows from the fact 
that the cutoff is larger than the VEV of S: S < M. A Gaussian perturbation spectrum, 
as indicated by observations, requires SS <C S, and therefore (5 < 1. For indirect decay 
5r = 3(s, whereas for the fluctuating mass case 5r = 3£s/2. 

The results for a non- negligible Q$ are shown in Fig. 2. The solid line corresponds to Q$ 
negligible small, when r( s <C (1 — r)£ 7 , we obtain ( = — 5 r /6 as before. The dotted line is 
for ( s = 10 2 and fi 5 w 10~ 5 , 2 x 1(T 5 , 5 x 10~ 4 at iV = 0, 4, 7 respectively. Such large 
curvature perturbations of S are possible for direct decay. The dashed line shows the plot 
for indirect decay, with £5 = 1/3 and Q$ ~ 5 x 10~ 3 , 0.1, 10 -2 at iV = 0, 4, 7. In both cases 
the final curvature is well approximated by Eq. ()4.12j) . with £ 7 = — 5r/Q- The contribution 
of (s to the total curvature rises above the one percent level for (^sCs)f(-y ~ 10~ 2 . For 
indirect decay and the fluctuating mass case, this is only possible for Q$ ~ 10 -2 . 



V. CONCLUSIONS 



In this paper we presented a gauge invariant formalism to study the evolution of curvature 
perturbations during decay of the inflaton or some other field that dominates the energy 
density. We specialized to the case where the perturbations arise from a spatially fluctuating 
decay rate. 

The system we considered consists of the inflaton <f), the radiation bath, and a light 
fluctuating field S which is responsible for the fluctuations in the decay rate of <fi. During 
inflaton decay there is a energy transfer from the inflaton field to the radiation bath. Due 
to the fluctuating decay rate this energy transfer has an intrinsic non-adiabatic component, 
which feeds the curvature perturbation of the radiation bath. This can be understood 
intuitively from the fact that fluctuations in the inflaton decay rate leads to fluctuations in 
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the reheat temperature of the Universe, and therefore to energy fluctuations in the thermal 
bath C 7 oc <5p 7 /p 7 oc ST/T. Since after inflaton decay the universe is radiation dominated, 
the total curvature perturbation is ( » £ 7 . 

We studied the evolution of the background fields and the perturbations numerically. We 
found that the total curvature perturbation is given by ( = — (1/6)8? and ip = (l/9)<5r- 
Here we have assumed that initially and Qs are negligible small. Departure from ( = 
— (l/6)5r, above the one percent level, arises for larger values of <^,(0)/£ 7 > 10~ 2 and/or 
fisCs/Cj ~ 10~ 2 . It is interesting to note that for indirect decay (s can be large: (s 3> <5r if 
the VEV of the flat direction field is much less than the cutoff scale S -C M. 

VI. NOTE 

While we were finishing our work we noticed a related contribution by Matarrese and 
Riotto, see jis] ]. In this paper, the gauge invariant curvature perturbations for 0,7 are 
presented. The authors have analytically obtained ip = (4/45)<5r/r, which in our case is 
proven numerically to be ip = (l/9)<5r/T. 
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